A companion to a quasar at redshift 4.7 
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t — There is a growing consensus that the emergence of quasars at high redshifts is related 
to the onset of galaxy formation 1 , suggesting that the detection of concentrations of gas 
accompanying such quasars should provide clues about the early history of galaxies. Quasar 
t> companions have been recently identified at redshifts up to z « 3 (refs 2-4). Here we report 
■ observations of Lyman-a emission (a tracer of ionised hydrogen) from the companion to 
a quasar at z=4.702, corresponding to a time when the Universe was less than ten per 
cent of its present age. We argue that most of the emission arises in a gaseous nebula 
• that has been photoionised by the quasar, but an additional component of continuum light 
Q\ perhaps quasar light scattered from dust in the companion body, or emission from young 
"*^h \ stars within the nebula appears necessary to explain the observations. These observations 
Q_i' may be indicative of the first stages in the assembly of galaxy-sized structures. 
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| The bright quasar BR 1202-0725 was discovered in the Automatic Plate Measuring optical survey for 
' quasars with z > 4 (ref. 5, see also ref. 6). The emission redshift is uncertain because Lya absorption 
produced by a foreground metal line system is present at a redshift z = 4.687 (ref. 7) that is much higher 
^ (by ~5000 km s" 1 ) than the redshift inferred from intermediate resolution spectroscopy 7,8 for the C iv 
$— i ' emission line of the quasar. 



We used the integral field spectrograph TIGER 9 (Traitement Integral des Galaxies par PEtude de leurs 
Raies) mounted on the 3.60 m Canada-France-Hawai Telescope during a windy run on 4-7 February 1995. 
The instrument samples the field of view through an array of micro-lenses, each providing a spectrum. 
Because of the windy weather during the observations, only two exposures were useful, corresponding to 
2.3 hours of integration. The seeing was 0.7" FWHM (Full Width at Half Maximum) and the instrumental 
configuration led to a spatial sampling of 0.61" in a circular field of view of 10" diameter centred on the 
quasar and a spectral sampling of 3.4 A. The spectral resolution, measured as the FWHM of an arc line, 
was 10 A, and the wavelength range was 5850-7100 A, corresponding to a Lya redshift range of 3.81-4.84. 
The data were reduced in a standard way. 10 ' 11 The la noise in the spectra obtained through one lens 
is 3xl0 -19 erg cm -2 A -1 after sky subtraction. We have then reconstructed 60 A wide narrow-band 
images. 

The line of sight towards BR 1202-0725 contains a known damped Lya absorption system with an 
H I column density iV(H i) ~ 4xl0 20 cm -2 at z = 4.383 (refs 7,12). Such systems are generally thought 
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to arise in proto-galactic discs 13 and so indicate potentially intensive star formation regions. However, 
we do not detect any emission down to a 3a limit of 3. 1CT 18 erg cm~ 2 A _1 arcsec -2 at z = 4.383 in a 
region of 41x41hj^ kpc 2 (where ^75 is the Hubble constant in units of 75 km s _1 Mpc -1 and q Q = 0.5) 
centred on the quasar. The limit reached, consistent with previous studies, confirms that damped Lya 
systems are not strong Lya emitters 14 ' 15 , either because of their small but significant dust content 16 or 
their low duty-cycle of massive star formation 17 . 

In contrast, the image centred at 6940 A (see Fig. 1) exhibits a significant (flux greater than 3d) 
North- West extension (region B) with the maximum emission located 2" away from the quasar (or 
9/175 kpc) that coincides spatially with an emission excess previously seen in broad-band VrI ls K 19 and 
narrow-band filter images 20 . The contribution of the quasar has been carefully subtracted, and we have 
coadded the spectra of all lenses in which the extension was detected (Fig. 2). There is a strong line (7 a) 
centred at 6932 A. The total flux, integrated over the line, is 2.5±0.6xl0 -16 erg s^ 1 cm~ 2 . Identifying 
the line with Lya gives a redshift z = 4.702. The emission is spatially resolved and has a maximum 
north-east/south- west extension of 1.5" after correction from seeing effects. 

As we detect only one line, other potential identifications could be Ha at z = 0.056, H/3 at z = 0.43 or 
[On]A3727 at z = 0.86. The probability of gravitational amplification is greatest in the latter two cases 
but population synthesis models indicate that the broad-band observations are inconsistent with these 
redshifts. In particular, the r — K colour characterises a steep spectrum that is difficult to reconcile with 
the lower limit on the B magnitude. Strong absorption systems are present at z = 1.7544 and 2.4421 
(refs. 7,12) in the quasar spectrum, but these redshifts are ruled out since no emission line is expected 
at 6932 A. 

The detected line coincides in wavelength with the peak of the quasar emission. This is a strong 
indication that the line is, in fact, Lya emitted by a quasar companion. Moreover, the blue wing of the 
line is sharp and coincides almost perfectly with the position of the first Lya absorption line seen in the 
quasar spectrum (see Fig. 3). This coincidence can be accounted for if the same absorber extends from 
the quasar to region B, implying a projected size of at least 9h^ kpc. 

The Lya flux is 1.7 times larger than the value obtained from previous narrow-band filter imaging 20 . 
A thorough check of the photometric quality of our observations using several standard stars yields an 
uncertainty of less than 20%. In addition, the flux determined for the quasar is consistent with published 
data 8 . The line luminosity is 2.1 xlO 43 /^ 2 erg s _1 and contributes 70% of the r-band emission. 

The observed flux is similar to that observed for quasar companions at lower redshifts (refs 2,3,21) 
and much larger than the upper limits set by surveys for primeval galaxies in blank fields 22,23 . This 
suggests that the presence of the quasar induces a peculiar interaction with the companion, either in 
the form of enhanced star formation or because the ionising flux from the quasar contributes to the 
excitation of the gas. In fact, the ionising flux of the quasar, estimated from the observed R magnitude 
and assuming a standard spectral index a = —1.5, is about 9.3xl0 56 photons per second. Under 
case B recombination, the observed Lya emission of the companion then implies that the fraction of the 
quasar luminosity intercepted by the gas is about 7xl0~ 4 . If we approximate the cloud as a sphere of 
diameter 6 kpc, corresponding to a covering factor of 3xl0~ 2 , the filling factor of the gas within the 
cloud is found to be over 0.02. The true value could be even larger if some of the Lya photons did not 
escape from the nebula. 

The interpretation of the extended emission near the quasar is constrained by the wavelength and 
equivalent width of the emission line, its contribution to the r band flux, the upper limit on the S-band 
emission, and the four VrIK broad-band fluxes. These constraints taken simultaneously rule out the 
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possibility that the emission line is not Lya at the redshift of the quasar. We consider hereafter two 
extreme models: a starburst galaxy and a gaseous nebulosity photoionised by the quasar. Examples of 
both situations exist at lower redshifts (z ~ 3; refs. 2,3,4). Virtually any model galaxy in which the 
spectrum is dominated by young stars is a reasonable match to the data. Fig. 4a shows a typical fit 
obtained for a 1.1 x 10 7 yr old model galaxy with metallicity 0.1Z@ and a constant star formation rate of 
13 M & yr" 1 . In this case, the onset of star formation would have occurred at Zf w 4.79 (again for ^75 = 1 
and go = 0.5). However, in starburst galaxy models the fraction of the r-band flux contributed by Lya 
emission is always less than 50 %, i.e., much smaller than observed. Therefore such models alone cannot 
explain the data. 

Alternative models of a pure photoionised nebula also present some shortcomings. If the / and K 
luminosities are assumed to be produced by reprocessing of the quasar's ionising radiation alone (Fig. 46 
shows such a model in which the reflecting cloud has a metallicity of O.1Z ), the predicted Lya emission 
is so strong that acceptable fits of the r band flux can be obtained only after assuming that a substantial 
fraction of Lya photons are removed from the line of sight, because, for example, of orientation effects 
or absorption by dust. 16 Arguing for a small filling factor of the cloud (see above) would not affect this 
conclusion because this would imply similar changes in both the Lya and continuum fluxes. The [O II] 
A3727 flux predicted for the model in Figure 46, 5x 10~ 18 erg cm~ 2 s _1 , is consistent with previous limits. 6 
The strong predicted N v emission could have been marginally detected (assuming that the cloud has 
zero metallicity would remove not only N V but also C IV emission, hence worsening the fit of the I-band 
flux). More important, in addition, is that the model cannot reproduce the observed V magnitude. 
Hence, although the Lya emission is most probably a consequence of the excitation by the ionising flux 
from the quasar, we conclude that an additional source of continuum is likely to be present. 

A possibility is that part of the continuum light from the quasar is scattered in our direction, poten- 
tially by dust. In fact, a large amount of dust has been detected in the quasar by continuum emission at 
millimetric and sub-millimetric wavelengths. 24 ' 25 New observations suggest that the emission is elongated 
towards the companion (A. Omont et al., mss in preparation), also indicating a large concentration of 
gas in that direction. We note that, by analogy with radio- galaxies, the ultraviolet continuum light of 
the companion could result from both scattered light and stellar emission. 26 The width of the Lya line 
(more than 1000 km s _1 , see Fig. 2) indicates strong disturbances of the kinematics of the gas. This all 
together suggests that we are witnessing the first stages in the build-up of a galaxy-sized structure. 
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Figure 1: Isophotes of the emission detected in the wavelength range 6910-6970 A plotted over the 
r-band image 18 . A Gaussian spatial filter was applied to the merged data cube, leading to a final spatial 
resolution of 1.2" (measured as the FWHM of the quasar image in the reconstructed images). The 
emission detected NW of the quasar peaks at 2" from the quasar. Contours of Lya surface brightness 
are 0.5, 1.2, 1.9, 2.6, 3.3, 4.0 and 4.7 10~ 16 erg s" 1 cm" 2 arcsec" 2 . 
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Figure 2: Spectra of the quasar (upper panel) and region B (lower panel). The position of a z =4.7 Lya 
line is indicated by a vertical dashed line. The quasar spectrum has been obtained by adding the spectra 
in all lenses in which its intensity is larger than three times the rms noise in the reconstructed broad band 
image. The quasar spectrum has been subtracted from the spectra of the lenses sampling region B. As 
the companion is detected only in the Lya line, the continuum map is an image of the unresolved quasar 
and can thus be used as the weight map to be applied to remove the quasar contribution everywhere in 
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Figure 3: Model fits of the emission line detected in region B (dashed lines) overplotted on the spectrum 
(solid line), (a) The model includes a single Gaussian profile absorbed in the blue by an intervening 
absorber. The emission profile peaks at A = 6931.5 A and has FWHM = 1245 km s _1 . The absorption 
line is centred at 6914.5 A and has FWHM = 6 A. The reduced x 2 is 1-04. (b) The model includes two 
Gaussian profiles centred at 6930.7 and 6951.4 A (separation 890 km s -1 ) and with FWHM = 900 km s^ 1 
and 410 km s _1 respectively. The reduced \ 2 is 0-96. 
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Figure 4: Model fits to the rest-frame emission properties of region B. Vertical error bars represent 
uncertainties in the broad-band observations, and horizontal error bars the filter bandwidths. Open 
circles represent filter-averaged model fluxes. Also indicated in each panel is the contribution of the 
Lya line to the integrated r-band flux, (a) Spectrum of a 1.1 x 10 7 yr old galaxy with metallicity 
O.1Z , a Salpeter initial mass function with lower and upper cutoffs 0.1 to 100 M Q , respectively, and 
a constant star formation rate of 13M Q yr _1 , as predicted by recent population synthesis models (G. 
Bruzual & S.C., mss in preparation). Nebular emission was computed 27 assuming case B recombination 
and gas parameters typical of giant Hli regions (n e = 100cm -3 , T e = 10 4 K, n(He n)/n(H)=0.1, and 
n(He m)/n(H) =0.01). The dotted line shows the resulting nebular continuum and the solid line the 
sum of the stellar and nebular continua with superimposed Lya and [O n] emission lines. Blanketing of 
the radiation blueward of Lya and photoelectric absorption of the (nebular) radiation blueward of the 
Lyman limit by foreground Hi absorbers at redshifts z < 4.687 have been included using approximate 
analytic formulae. 28 The dashed line shows the total (stars plus nebular) spectrum in the absence of Hi 
absorption, (b) Spectrum reflected by a cloud with metallicity 0.1Z Q photoionised by the quasar according 
to standard photoionisation models. 27 The ionising spectrum of the quasar has been approximated by a 
power law f v oc u -1 ' 5 . The electron density, Lyman limit optical depth, and ionisation parameter of the 
model are n e = 100 cm" 3 , tl = 100, and log U = —1, respectively, but the quality of the fit depends only 
weakly on these parameters. To improve the fit, Lya emission has been reduced arbitrarily by 90%. As 



